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I. INTRODUCTION 

The discussion of periodic motion in one dimension is important in most introductory courses on classical mechanics. 
Several problems can be solved exactly, but in most cases one has to resort to approximate solutions. Simple but 
sufHciently accurate approximate solutions for such problems are very important in understanding many features 
of classical mechanics. In addition to it, in some cases one is simply satisfied with accurate numerical results, and 
expressions suitable for computation are most welcome. 

The purpose of this paper is the discussion of the exact and approximate calculation of the period of a particle that 
moves in one dimension under the effect of an anharmonic potential. 



II. PERIODIC MOTION IN ONE DIMENSION 

Consider a particle of mass m moving in one dimension under a potential-energy function V{x). Without loss of 
generality we assume that V{x) has a minimum at x — 0; more precisely, we assume that V{Q) — 0, V [Q) — 0, and 
V"{Q) > 0, where the prime indicates differentiation with respect to x. Following the standard notation in classical 
mechanics, we use a dot to indicate differentiation with respect to time, for example: v = x. 

^From the equation of motion 

mx = -V'ix) (1) 

we easily obtain an integral of the motion 

E='^ + V{x) (2) 

which is the total energy. The motion of the particle is restricted to the interval x- < x < x+, where the turning 
points x± satisfy V{x±) = E; that is to say, i; = at those points. 
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It is well-known that the period of the motion is given by 

(j>dt = V2^ , "'^ (3) 

J 1_ ^E-V{x) 

from which we obtain the frequency = 27r/T. 

We can simplify the equations of motion by the introduction of a dimensionless time r = tjjf)t, where wq is an 
arbitrary frequency. If we define 

S ^ A, f^(-) = ^ (4) 
mujQ mujQ 

then we obtain the equations of motion of a particle of unit mass; for example: 

£^^ + Uix) (5) 

and 

^0 ^£ - U{x) 

It is worth noticing that equation 101 is not dimensionless because £ and U{x) have units of length squared. In 
order to get a truly dimensionless equation we should define a dimensionless coordinate q — x/L, where L has units 
of length. Thus E/ {itiujqL'^) and V{Lq)/ {mujQL^^ are the dimensionless counterparts of the total and potential 
energies, respectively. In this paper, however, we have opted for equations that are similar to those often found in 
current literature. 



For example, if 



2 4 
V2X V4X 



we may choose ojo = \fv2lm so that 



where A — W4 /w2 • 



U{X) = y + — , (8) 



III. THE MAIN INTEGRAL 



It follows from the discussion above that the period is proportional to an integral of the form 

dx 

1= / -f^, (9) 

were Q{x) exhibits simple zeros at x^ and xj^ and is positive definite for all X- < x < xj^. That is to say, we can 
write 

Q{x) = - x)(x - x^)R(x) (10) 

where Ti(x) > for all x_ < x < x+. 
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The reason for rewriting our problem in this somewhat abstract way is that the integral ® applies to problems 
other than the period of a motion in one dimension. We will mention some examples later on. 

In order to develop suitable exact and approximate expressions for the integral ^ we define the reference function 

^2 

Qo(a;) = — (a;+ - x) (a; - (11) 

that satisfies the appropriate boundary conditions at the turning points. It is clear that Qo{x) is the function that 
would appear in the treatment of a harmonic oscillator. Then we rewrite (O as 

/- '^'^ (12) 



where 



The change of variables 



makes the integral l|12|) much simpler: 



Q(x)-Qo(x) 2R{x) - u? 



X = -i— + — cose (14) 



Jo VlTA 

This equation leads to an exact expression for the period, which in most cases one has to calculate numerically. In 
order to derive simple analytical formulas we expand 



where | | = al/[bl{a — b)l] is a combinatorial number. Notice that this series converges for all x such that |A| < 1. 
We thus obtain a series for the integral (|15|l : 



r'^' I / ^^do. (17) 



In this way we can derive approximate expressions for the integral (|15() by means of the partial sums: 

N 

(18) 

IV. THE DUFFING OSCILLATOR 

The potential-energy function © gives rise to the Duffing oscillator. Since it is parity invariant (U{—x) = U{x)) 
then x-i- = —X- = ^4 is the amplitude of the oscillations. According to the general discussion of the preceding section, 
it follows from 



Q{x) = £ - U{x) = {A^ - x^) 



(19) 
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that 



and 



Rix) = i + ^ (A^ + x^) (20) 



A = 5 ^ 21 

where we have substituted x — A cos 0. We conclude that the period depends on the dimensionless parameter p = XA^ 
that is the ratio of v^A"^ and V2A^ both having units of energy. 

If we choose lo = yT~Fp then we obtain an ah'cady known suitable compact expression for the integral^ 

V2 r d0 



'^=7^-^- (22) 



This equation yields the series 




^=^^1. e (23) 



that converges for all |^| < 1; that is to say, for all p > —2/3 or p < —2. 

When A < the potential exhibits two barriers of height l/(— 4A) at a; = ±l/\/— A and therefore the amplitude of the 
periodic motion cannot be greater than A^ = l/\/— A. In other words, there is periodic motion if p > = AA| = — 1. 
The series H23|l does not converge for —1 < p < —2/3 and the analytical expressions that we may derive from it will 
not be valid for all the values of the energy that give rise to periodic motion. Can we improve this approach?. The 
answer is "yes " as we will see below. 

Let Rm and Rm be the maximum and minimum values of R{x) in the interval and Am and Am the 

corresponding values of A(x). Since R{x) is positive definite we know that Rm > R{x) > i?,„ > 0. If we choose the 
value of the adjustable parameter lo so that Am ~ ~A„i we obtain 

u;,2 = Rm + Rm>0 (24) 

and 

2R{x) - Rm - Rm 

^b[x) = — . (25) 

Km + n-m 

The subscript b indicates that this particular value of uj "balances" the maximum and minimum values of A (a;). 
Notice that |A(x)| < 1 for all x^ <x < x+ because Aj\/ = {Rm — Rm) / {Rm + Rm) < 1- 

For the particular case of the Duffing oscillator we have R„i — R{0) — 1/2 + p/A and Rm — R{±A) — 1/2 + p/2 so 
that 

^f, = — : — (26) 



and 



Ah(e) = j^cos(20). (27) 
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Thus the integral becomes 

2^/2 r de 



V^+^Jo ^l + ecos(20)' 4 + 3p 

that gives rise to the series 



(28) 



'""'.ti-^y -'n -"'\e'. (29) 



and 



which converges for all |^| < 1; that is to say, for all p > — 1 or p < — 2. In this way we may obtain simple analytical 
expressions for the period valid for all values of the energy consistent with periodic motion. 
According to equation © the period is given by 

T = —I (30) 

UJq 

and equation 1)29(1 enables us to derive simple analytical approximate expressions for it. For concreteness and simplicity 
we choose luq ^ 1 in what follows. For example, the first two approximations are 

r(o) = -i^ (31) 

^ ^(l47p^ + 384p + 256) 
4(4 + 3p)'/' 

These expressions are expected to be accurate for small values of p, and in fact they give the exact result for p = 0. 
However, they are also accurate for extremely great values of p. Notice that 

lim = 4 / , w 7.4162987 (33) 

^ Jo ^3 + cos(20) 

A straightforward calculation shows that 

lim y^T(°'> - ^ « 7.26 (34a) 

p^oo ^3 

lim VPT^'^ - ~ 7.406. (34b) 

p^oo 36 

We conclude that such simple analytical expressions are sufficiently accurate for most purposes and that one easily 
improves them by straightforward addition of more terms of the series ()29f) . 



V. QUADRATIC-CUBIC OSCILLATOR 

Parity-invariant oscillators exhibit symmetric turning points; if the potential is nonsymmetric so are the turning 
points. The simplest example is 

nx) = fx^ + |x^ (35) 



If we again choose ujq = ^/v^jm then we obtain 



U{x)^i- + \x\\^''±. (36) 
2 3 V2 
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The potential-energy function U{x) shows a barrier of height U{xm) = 1/ (6A^) at xm = ^ and the turning 
points satisfy x+ > > x^. 
If we write 

Q{x) = £ - U{x) = (x - x-){x+ - x){bo + bix) (37) 

then we obtain 

x+x^ A X++X- 

^ .2 I ^ ^ — rzTT ' ^1 = - = -T^TZT—: — -3T (38) 



2{x^_^ + x+x^ + x±) 3 2{x^ + x+x^ + xi) 



and 



3 

2 + x+x^ + x\ 

Since U{—x, —A) = C/(a:^, A) we consider only the case A > without loss of generality; therefore X- + x+ < because 
+ x+X- + x'L > 0. Taking into account that 6o > and 6i > we conclude that Rm and Rm take place at the 
turning points; therefore, 

n , . , . _l_ I '^X -4- X _ I ijj 

Lol = + i?(x_) = + ^ + (40) 

2 (a;!j_ + a;+a;_ + x_) 

Is real for all values of £ below the barrier?. In order to answer this question notice that the third root X3 of Q{x) 
is smaller than x_ and is given by 

hr, x+x^ 

X3 = -7^ = ^<0. (41) 

Therefore 



x'^ + 4x+a;_ + x'^ 



— {x^ + x^)'^ + 2x^x^ 

= {x++x-){x++x--2x3) <0 (42) 

because a;_ — 2:3 > and X4- — X3 > 0. 

Finally, after the change of variables (|14|) the function A{9) takes a particularly simple form: 

A,(ff)^gcos^, , 2 (43) 

where x'^ — xt < because < a::+ < — 
The resulting integral 

/=^r^^^ (44) 

Jo Vl + Ccos6' 



gives rise to the series 



■ n 



3 J \ 2i 

which is similar to the one derived above for the Duffing oscillator and converges for all |^| < 1 
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When £ = U{xm) then X3 = a;_ (remember that A > 0) and ^ = 1. We appreciate that the series (|45(l converges 
for all values of the energy for which there is periodic motion. 

Equation H44|l gives us a simple and exact expression for the period of the anharmonic oscillator that requires 
numerical integration to obtain results for a given set of potential parameters. On the other hand, equation H45() 
provides approximate analytical expressions that one makes as accurate as desired by simply adding a sufhciently 
large number of terms. The choice of one or another depends on the particular application. 

Following a different procedure Apostol^^ derived the exact expression for the period of the quadratic-cubic oscillator 

, ^ r^' = (46) 

uja^/x+ - X3 Jo ^1 _ fc2 gin2 a X+-X3 

The expansion of this equation in powers of fc^ is also convergent for all values of the energy consistent with periodic 
motion because < 1. 




VI. CONCLUSIONS 



In this paper we present a straightforward systematic procedure for constructing exact and approximate expressions 
for the period of anharmonic oscillators. The recipe is simple: first, we factor the function Q{x) and obtain the turning 
points and the function R{x) as in equation H1U|) . Second, we obtain the maximum and minimum values of R{x) in 
the interval between the turning points which determine the optimum value of uj. Thus we are left with an exact 
expression for the period that we may use in numerical applications. In addition to it, we may expand this exact 
expression in a Taylor series in order to obtain partial sums that become analytical expressions for the period of 
increasing accuracy. These partial sums converge to the exact result for all values of the energy that give rise to 
periodic motion. 

The method proposed in this paper is not restricted to the period of anharmonic oscillators with polynomial 
potentials. We may, for example, expand a given arbitrary potential U{x) about its minimum to any desired degree 
and then apply the approach developed above. Moreover, some other problems have been expressed in terms of 
integrals of the form such as, for example, the deflection of light by a massive body or the precession of a planet 
orbiting around a star—. Recently, we have already applied a variant of present approach to such problems^. 

There is a wide range of interesting applications for present method and for that reason we believe that it is suitable 
for teaching in advanced undergraduate courses on classical mechanics. 
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